Spin-entangled currents created by a triple quantum dot 
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We propose a simple setup of three coupled quantum dots in the Coulomb blockade regime as 
a source for spatially separated currents of spin-entangled electrons. The entanglement originates 
from the singlet ground state of a quantum dot with an even number of electrons. To preserve 
the entanglement of the electron pair during its extraction to the drain leads, the electrons are 
transported through secondary dots. This prevents one-electron transport by energy mismatch, 
while joint transport is resonantly enhanced by conservation of the total two-electron energy. 
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The creation of entangled particles is a crucial prob- 
lem as entanglement is a prerequisite for quantum com- 
putation and communication Jj]]. While manipulations of 
entangled photons demonstrating various quantum infor- 
mation processing schemes have been very successful 
similar achievements are still missing for massive parti- 
cles such as electrons. Hence there has been a number 
of theoretical proposals for a solid-state entangler - a de- 
vice creating two entangled particles and allowing their 
separation and extraction into two distinct channels for 
further processing. 

Recent proposals involved the extraction of entangled 
Cooper pairs of a superconductor in contact with quan- 
tum dots ||, normal or ferromagnetic conductors (|, |j, 
and carbon nanotubes || [jj. In another scheme, the en- 
tanglement arises from interference effects in a quantum 
dot in the cotunneling regime and requires special non- 
degenerate leads of narrow energy width |J. A generic 
entangler based on interferometry and which-way detec- 
tion was proposed in Ref. ||. In this article, we propose 
an entangler based on a triple quantum dot setup. The 
entanglement originates from the singlet state of a pair of 
electrons in one quantum dot, while its transport relies 
on energy filtering by secondary dots. Our proposal is 
based on existing technology (To) and on realistic param- 
eter values as typically found in transport experiments 
with quantum dots. 

Setup. Fig. |j] describes the proposed entangler. It is 
composed of three coupled lateral quantum dots (Dc, Dl 
and Dr) in the Coulomb blockade regime, each of them 
coupled to a Fermi liquid lead lc, II and Ir. When two 
excess electrons are present in Dc, we can assume Jll[] 
that their ground state is the spin-singlet state, which is 
the (anti) symmetric superposition of their (spin) wave- 
functions. The aim of the entangler is to extract the 
singlet from Dc, by transporting one electron into the 
neighboring dot Dl and the other one into Dr, and fi- 
nally transport them into the drain leads II and Ir with- 
out loss of entanglement. This creates two currents of 
pairwise spin-entangled electrons that are spatially sep- 
arated. 

Applying two bias voltages fic — fJ-L and fic — al- 
lows the transport of electrons from the source lead lc 




Figure 1: (a) Setup of the triple quantum dot entangler. Three 
leads li (i = C, L, R) at chemical potential m are coupled to three 
quantum dots D; in the Coulomb blockade regime. Each dot con- 
tains an even number of electrons, and can only accept 0, 1 (Dl 
and Dr) and 2 (Dc) excess electrons. A spin-singlet is formed 
in Dq when two electrons tunnel incoherently, each with a rate 
a, from the source lead lc into Dc- Each of the electrons can 
subsequently tunnel coherently to Dj, and Dr with tunneling am- 
plitudes Tl and Tr. Finally, the electrons tunnel out (with rates 
7_L and 7_r) to the drain leads II and Ir, creating two currents of 
entangled electrons, (b) Energy level diagram, for each electron. 
Non-entangled currents, which arise from one-electron transport, 
are suppressed by the energy differences £l,r — ec ±U. The joint 
transport of both electrons is favored by conservation of the total 
two-electron energy: el + er ~ 2ec (resonance condition). 

to both drain leads II and Ir, via the three quantum 
dots Di,i = C, L, R. To preserve the entanglement of 
the electrons until they are in the drain leads, one must 
avoid the individual transport of one electron, as this 
would allow the arrival of a new electron in Dc which 
could destroy the existing entanglement by forming a 
new singlet with the remaining electron. To suppress 
one-electron transport, we arrange the dots so that there 
is a large difference between the energy levels of Dl and 
Dc compensated by the energy difference between Dr 
and Dc- This way the joint transport of both electrons 
to each neighboring dot conserves energy and is therefore 
enhanced by resonance, while the off-resonant transport 
of one electron is suppressed by the energy difference. 

The number of electrons participating in the transport 
is controlled via Coulomb blockade j[2), where N excess 
electrons in dot Di create a large electrostatic Coulomb 
charging energy Ui(N). The energy of the N th electron 
is then Ei(N) = U l (N)~U,(N - l) + e l (N), where e^N) 
is the lowest single-particle energy available for the A^*' 1 
electron. We consider the electrons to be independent 
and neglect further effects such as interdot charging en- 
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ergy or exchange Coulomb interaction. Assuming a shell 
filling of each dot [|l2| , we disregard all but the excess elec- 
trons in each dot. The ground state in D c with 2 excess 
electrons is the spin singlet f-l — IT, where both elec- 
trons have the same orbital energy ec — ec(l) = e c(2). 
One excess electron in Dl or Dr cannot form a singlet 
with one of the electrons already present as these are 
already all paired up in singlets. We assume the gate 
voltages of each dot so that Uj(0) — Uj(l),j = L,R and 
Uc(0) = Uc(2) — 0, which gives a negative charging en- 
ergy for one electron: £7c(l) = — U. We define the zero 
energy as the total energy of the three empty dots. 

The energy levels in the dots and the chemical poten- 
tials /Ltj in the leads are assumed to be tuned such that 
only zero or one excess electron in Dj,j = L,R, and 
zero, one, or two electrons in Dc are allowed. It is cru- 
cial that only the ground states of the electronic levels 
in the dots participate in the transport. In particular, 
the triplet states in Dc should not be accessible to in- 
coming electrons. To avoid resonance with excited lev- 
els, the energy level spacings in the dots must be larger 
than the Coulomb charging energies: Ae^ > XJ{. Excited 
states with energy E* could participate in the transport 
through cotunneling events p2[ , where one describes the 
transport of one electron from lc to lj via intermediate 
"virtual states" in the dots by second-order processes. We 
shall neglect such events, as they are suppressed by fac- 
tors of the order of a/ (E* — /Lt«) ~ a/(Ui — fa) <C 1. 

We need low temperatures T so that thermal fluc- 
tuations cannot allow three electrons in Dc or popu- 
late excited levels, which could also create a current in 
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the reverse direction (lj 



Dj, D c lc). Taking 



k^,T <C \fa — Ei(0, 1, 2, 3)1, we can neglect temperature ef- 
fects and set T = for simplicity. The incoherent tunnel- 
ing rates are given by a = 2ir\tc\ 2 vc and jj = 2-K\tj\ 2 Vj, 
where ij is the tunneling matrix element connecting U 
and Di, and i/< is the density of states of the lead. 

The quantum states of the entangler are given by 
combining the different numbers of electrons allowed in 
each dot. describes the situation where all the dots 
are empty; L, R or C corresponds to one electron in 
Dl, Dr or Dc, while CC denotes the singlet state 
created by two electrons in Dc- Thus, the 8 states, 
shown in Fig. ^| with their transitions, for the basis set 
B = {0,C,CC,LC,CR,LR,L,R}. This description in 
terms of the individual levels of each isolated dot requires 
that the tunneling matrix elements Tl and Tr (consid- 
ered to be real) connecting the dots are small and do not 
mix the levels in different dots, i.e. Tl,Tr <C Ae*. 

The coherent evolution is described by a Hamiltonian 
matrix Hk t k'>k,k' <E B. The diagonal elements contain 
the energies E — 0, E c — e c - U,E C c = 2e c ,E L c = 
eL+ec - U, E C r = cr + cc — U, E LR = e L + e«, E L = 
£l,Er = €r, while the off-diagonal elements describe 
the coherent oscillations of one electron between Dc and 
Dl or Dr: Hc,l = Hcr,lr — Tl,Hc,r = Hlc.lr — 
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Figure 2: The entangler states {0, C, CC, LC, CR, LR, L, R} and 
their transitions. Double arrows indicate coherent tunneling (os- 
cillations) of one electron between two dots with overlap matrix 
element Tl and Tr (in the shaded area). Single arrows indicate 
incoherent tunneling from (a) or to (jl and 7_r) the leads. The 
dashed lines indicate the three types of transitions that must be 
avoided to ensure the joint transport of the singlet pair CC to the 
leads (see text; for clarity we do not show the transitions obtained 
by replacing L by R.) 

Tr, H cc ,lc = T L s/2 and H cc ,cr = T R V2. The \/2 
factor comes from the identical orbital states in CC. 

We describe the incoherent transport as sequential 
tunneling (lowest order in 0,7^,7^) in terms of the 
master equation for the (reduced) density matrix p 
of the entangler. The diagonal elements, pk,k S 
B, are the occupation probabilities of the state k, 
with normalization J2k Pk = !• The off-diagonal ele- 
ments pk.k' contain the coherent superposition of k and 
k! . We write the master equation [jl3| as dp /At = 
— i[H,p] — M, where M describes the incoherent trans- 
port from/to a lead connecting a state k to another 
state k' with a rate W(k',k) 6 {a, 7i,7ij}. For the 
diagonal elements, this results in the population equa- 
tion: M k = ~PkE k 'W{k',k) + J2 k ,W{k,k')p k ,, while 
the off-diagonal elements are damped by the incoherent 
transitions out of each of the two corresponding states: 
M k>k , = -\ PKk , J2 k »^ k , k > W{k", k) + W(k", k'). Out of 
the 64 elements of p, only 24 are coupled to the rele- 
vant diagonal elements. Arranging them in a real vec- 
tor V , one can rewrite the master equation as a homo- 
geneous first-order differential equation given by a ma- 
trix A: dK/dt = AV . Its stationary solution is the 
eigenvector corresponding to the zero eigenvalue of A. 
It can be found symbolically by a mathematical soft- 
ware (MAPLE), and defines the stationary populations 
P k = p k of the different states k of the entangler. 

Results. Our aim is to show that the transport of the 
singlet state through the entangler is robust and domi- 
nates over the transport of uncorrelated electrons. We 
define the stationary entangled current (see the central 
sequence in Fig. |) as I E = e^ L (P LR + P L ) + ej R (P L R + 
Pr). The destruction of the entanglement of the two elec- 
trons can occur in three ways; see Fig. || (the following 
discussion also refers to the cases with L replaced by R) . 

(1) An electron tunnels out from Dr to Ir while the 
second electron is still in Dc- This creates a current 
z-IrPcr which might contain no entanglement, as the 
remaining electron can form a new singlet with a new 
electron coming from lc- This current is suppressed by 
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the energy mismatch cr — ec — U between the states CC 
and CR. The entangled current Ie is dominant if the 
entangler quality Q = rmn{P LR / Pgr, PlrI Plc} satis- 
fies Q > 1. 

(2) An electron tunnels from Dc to Dl and then out 
to I l before a second electron has tunneled into Dc 
and formed a singlet. This channel is suppressed by 
the energy difference e^ — ec + U between the states 
C and L. As the stationary solution contains no in- 
formation on the history of the electrons, we compare 
currents obtained in the two following cases, (i) We 
switch off the undesired channel by setting Tj, = Tr = 
between C, L and R, while keeping the tunneling be- 
tween CC,CR, LC , and LR. This defines the station- 
ary populations Pk- (ii) We keep the undesired chan- 
nel, while switching off the tunneling involving CC and 
LR. This defines the populations Pk, and creates a 
current e"/ L (PL + Plc) + ej R (PR + P C r) containing 
no entanglement. Defining a second entangler quality 
by Q = mm{P L /PL,PLR/PLC,PR/PR,PLR/Pc R }, the 
condition Q> 1 corresponds to the suppression of this 
one-electron channel. 

(3) After the joint transport of the two electrons into 
state LR and the tunneling of one electron into It, a new 
electron tunnels into Dc before the remaining electron in 
Dr has tunneled out to Ir. The new electron can then 
form a new singlet with the remaining electron, there- 
fore destroying the entanglement that existed with the 
electron which has moved to the lead II- To suppress 
this channel, we need a <C 7l,7_r- For simplicity, we set 
a = for transitions from L and R when calculating the 
probabilities Pk (but we keep a non-zero for Pk). 

The exact analytical expressions for Pk are extremely 
lenghty and cannot be written in a compact form. We 
introduce Sej = e 3 — ec, which explicitly removes the 
dependence on ec as only energy differences enter the 
time evolution of the density matrix p(t). Secondly, we 
consider a symmetric setup with 7 = — TR and 
To = Tl = Tr. Then, the symmetry between LR, L, 
and R yields P L = P R = P LR (see Fig. |). Thirdly, we 
consider the case €r = ec- The resonance condition for 
the joint transport Ecc = 2ec — Elr, = (l + £r trans- 
lates into SeL — 0, and the energy differences relevant for 
the suppression of one-electron transport are given by U 
and U ± SeL- From a qualitative analysis involving first 
and second-order perturbation calculations (a la Fermi 
Golden Rule), we find that the conditions for a domi- 
nant joint transport read a <C 7 <C T <C U,\U ± 5ei,\ 
|p^| . Hence we expand both the numerator and denom- 
inator of the expressions for Pk in the lowest non-trivial 
order in a, 7, To. We distinguish two cases: 
[I] non-resonant (\Sei,\ > To): 
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[II] resonant (Scl = 0): 
Pl=Pr = Plr 
Pgr = Plc - 
Pr = Pl 
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and P C R = P R a/j, P LC = Plcy/i- 

These two cases are sufficient for an analytical discus- 
sion of the entangler, as the approximate expressions (|l|)- 
@ reproduce accurately the exact results presented in 
Fig. |^. As shown in Fig. ||(a), the entangler quali- 
ties Q and Q reach a maximum around SeL = 0, which 
is due to a resonance in the coherent oscillations be- 
tween^ CC and LR. We define now the quantities Q\ 
and Q\ as approximations of the qualities Q and Q ob- 
tained with Eqs. ([j])-(||)- Qi and Q\, which grow as 
(Se^)^ 2 , give a correct estimate of the width Se^ of the 
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Figure 3: Quality of the entangler and output current, with 
the parameters a — 0.1,7 = l,To = 10, U = 1000 in /LteV. 
(a) Quality Q and Q, around the resonance at SeL = where 
the entangled current dominates. The width of the resonance 
defined by Q, Q > Qf n = 10 is \5e L \ < Jk~l = 2T / A /Qf in ~ 
6 /jeV (shown in gray), (b) Q and Q as a function of the 
tunneling matrix element To at resonance (SeL =0). In gray, 
the region where the quality of the entangler is Q, Q > Qn' n 
1 UO. (c) Average current in the left {IT) and in the right 
{Ir) leads, which takes into account both entangled and non- 
entangled currents. A symmetric current (J£ v = Iff) is the 
signature of the resonance and therefore of the desired regime 
where the entangled currents dominate (see the inset with 
a larger scale), (d) At resonance, the entangled current Ie 
saturates to ~ a = 4pA when T„ > {U 2 -ya/32) 1/4 ~ 5peV. 
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resonance around Scl — 0. Introducing the condition 
QuQi > Qf n , we get \Se L \ < Je~ L = 2T / \/Qf™- Note 
that such proportionality of the width to the tunnel- 
ing matrix element is also found in the Rabi formula 
for a two-level system. Similarly, we define Qu and 
Qn with Eqs. these approximate ratios accu- 

rately reproduce the height of the resonance peak at 
SeL = 0. Introducing the condition Qn,Qu > Qu m , 
we find 7\/ ( 2n I 7 8 < T o < u V 4a /7<2n m ; see Fig. |(b). 

Quantities which are experimentally accessible are the 
currents in the left and right leads, II — £Jl(Plc + 
Plr + Pl), and I R = e-f R (P C R + Plr + Pr)- Far from 
resonance (SeL > 2q), one cannot neglect the currents //, 
and I R coming from the single-electron transitions de- 
scribed by the probabilities Hence we consider in 
Fig. |^(c) the average of the currents in situations (i) and 
(ii): Tf = [I L + I L )/2 and 1% = [I R + I R )/2. Away 
from the two-electron resonance the current is asym- 
metric: Il v ^ 1%. The large peaks where J£ v > 1% 
are due to one-electron resonances: between CC and 
LC when SeL = U (right peak), and between C and 
L when SeL = —U (left peak). At the SeL — res- 
onance the current is symmetric as the electrons are 
transported together and simultaneously from the cen- 
tral dot to the leads [see inset in Fig. §|(c)]. Hence by 
varying e_L via the gate voltage of Dl until II = I R , 
one can locate the resonance, and therefore the regime 
where the entangler is most efficient Finally, Eq. 

(|) gives P LR -> a/4 7 for T 4 > [/ 2 7 a/32, which yields 
Ie —> W = a, as illustrated in Fig. ||(d). Note that 
in order to be able to reach I n 



within the window 

1\/Qn n / 8 < T o < lV a / 4 7<2n in > one also needs the 
condition U > Q^ in -/y / j/2a. 

Discussion. Setting the qualities of the entangler to 
Qf™ = 100, Qf in = 10 and 7 = 10a, we need, approxi- 
mately, 35a < T < C//60, and U > 2200a. The first con- 
dition is easily met as a and To can be varied via the volt- 
ages defining the barriers For the second condition, 
present-day experiments manipulate currents typically 
around 2pA 0, yielding a > O.l^cV^C/ > 0.3 meV. 
To get a finite width SeL — 6/xeV, we take U ~ lmeV, 
which is within reported values |jl^, [ll|. To obtain an 
even better quality (Qjjf 111 = 1000,7 = 100a), one needs 
to increase the ratio U J a. One possible issue is the single- 
particle energy spacing Ae^, which is usually smaller than 
or equal to the charging energy. However, as Aei oc l/L 2 
in a box of size L, while U oc l/L, one should be able to 
reach Aei ^> U by decreasing the dot size. Alternatively, 
one can use vertical quantum dots, which have large en- 
ergy level spacings Jl^| , or use one carbon nanotube with 
two bendings (which defines three regions behaving like 
quantum dots). 

A current I max ~ 10 pA corresponds to the delivery of 
an entangled pair every t p ~ 1/a ~ 16 ns, which is below 
reported spin decoherence time of 100 ns |0. The av- 



erage separation between two entangled electrons is ap- 
proximately t G — 1/U ~ 0.4 ps <§; t p , with a maximal 
separation of t m ~ 1/7 ~ 0.5 ns. This would allow noise 
measurements using a beam splitter, where an enhance- 
ment in the two-terminal noise is a signature of singlet 
states compared to the noise of non-correlated electrons 
fUjfl , Note that such an experiment requires electrons 
with the same orbital energy, which can be achieved if 
(l = (r — ec- One could also carry out a measurement 
of Bell's inequalities by measuring the spin of the elec- 
tron in each lead with the help of a spin filter based on 
a quantum dot Jl9| |. 
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